arXiv:1509.02774v3 [nlin.PS] 17 Jul 2016 
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While the existence of stable knotted and linked vortex lines has been established in many exper¬ 
imental and theoretical systems, their existence in oscillatory systems and systems with nonlocal 
coupling has remained elusive. Here, we present strong numerical evidence that stable knots and 
links such as trefoils and Hopf links do exist in simple, complex, and chaotic oscillatory systems if 
the coupling between the oscillators is neither too short ranged nor too long ranged. In this case, 
effective repulsive forces between vortex lines in knotted and linked structures stabilize curvature- 
driven shrinkage observed for single vortex rings. In contrast to real fluids and excitable media, the 
vortex lines correspond to scroll wave chimeras [synchronized scroll waves with spatially extended 
(tubelike) unsynchronized filaments], a prime example of spontaneous synchrony breaking in systems 
of identical oscillators. In the case of complex oscillatory systems, this leads to a novel topological 
superstructure combining knotted filaments and synchronization defect sheets. 

PACS numbers: 05.45.Xt, 89.75.Kd, 82.40.Bj, 02.10.Kn 


Introduction: In natural science, knots and linked 
structures have attracted attention in various branches 
as they are an essential part of many physical processes. 
This includes real fluids [T|, liquid crystals mm, Bose- 
Einstein condensates H E], electromagnetic fields and 
light [H |7], superconductors [8], proteins [9] as well as 
excitable media mm and bistable media [12]. Stable 
knots and their topological invariants are of particular 
interest for both theory and experiments as they play an 
important role in characterizing and controlling different 
systems and their dynamics [13]. This is especially true 
in excitable media, where linked and knotted filaments of 
phase singularities can be essential to understand the na¬ 
ture of scroll wave propagation processes [IU1III1II1II5], 
including nonlinear wave activity associated with ventric¬ 
ular fibrillation and sudden cardiac death [lailT]. While 
the wave propagation processes in excitable and nonlin¬ 
ear oscillatory systems are very similar [18], the existence 
of such stable knotted and linked filaments in oscillatory 
systems has remained elusive. For example, to the best 
of our knowledge no corresponding parameter regime has 
been identified in the compiex Ginsburg-Landau equa¬ 
tion (CGLE), which is the normai form of osciiiatory 
media ciose to the Hopf bifurcation [191 EQ]. This is 
deepiy unsatisfying as the coiiective behavior, sponta¬ 
neous synchronization and wave propagation in osciiia¬ 
tory media and coupied systems of noniinear osciiiators 
are topics of generai interest with appiications across dis- 
cipiines [Hill], inciuding the quantum regime [2211^ . 

In this paper, we show for the first time that (i) stabie 
knotted and iinked fiiaments do exist in osciiiatory sys- 
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terns, (ii) they do exist under non-iocai coupiing in the 
underiying dynamicai equations, and (iii) together with 
synchronization defect sheets they can form novei topo- 
iogicai superstructures. From the Kuramoto modei of 
simpie phase osciiiators and the CGLE to compiex and 
chaotic osciiiatory systems, we find in particuiar Hopf 
iinks and trefoiis that persist over hundreds of thousands 
of scroii wave rotations for a wide range of parameters. 
Due to the non-iocai coupiing, the fiiaments that make 
up the iong-iived knotted structures are no ionger sim¬ 
pie phase singuiarities as is typicai for scroii waves, but 
instead the fiiaments correspond to spatiaiiy extended re¬ 
gions in which the osciiiators are unsynchronized. This 
is despite the fact that aii osciiiators are identicai and 
uniformiy coupied. The coexistence of these unsynchro¬ 
nized iocai regions with synchronized regions — exhibit¬ 
ing traveiing waves in our specific case — is the haiimark 
of a chimera state [MH29] . Whiie singie ringiike chimera 
fiiaments shrink, knotted and iinked fiiaments generate 
an effective repuision that prevents shrinkage and stabi- 
iizes the pattern even in the presence of strong noise. We 
find that for coupiing that is too short ranged (inciuding 
iocai coupiing) or coupiing iag that is too smaii, the re¬ 
puision is too weak such that knotted structures coiiapse. 
This is despite the fact that phase twists aiong the fiia¬ 
ments are present, which have been hypothesized to have 
a stabiiizing effect by themseives m- If the coupiing be¬ 
tween osciiiators is too iong ranged and the coupiing iag 
is too iarge, straight chimera fiiaments become unstabie 
in a way reminiscent of negative iine tension [niEoHas]. 
This ieads to the breakup of the knotted structures as 
weii. 

Phase oscillators: As the simpiest paradigmatic 

modei of an osciiiatory system, we first focus on the Ku¬ 
ramoto modei EHlEllET]. In this modei, 6>(r,t) G [—7r,7r) 
denotes the state of an osciiiator at a spatiai point r and 
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time t. The evolution is governed by 

^(r, t) = cjo + Kuj{r, t). (1) 

Here, ujq is the natural frequency of the oscillators and 
K is the coupling strength. Note that we can set cjo = 0 
and rescale time Kt ^ t without loss of generality. Thus, 
uj is the instantaneous angular frequency obeying 

uj(Y^t)= f Go{r — r') sm[0{r\t) — — a]dr' (2) 

Jv 

where Go(i*) is a coupling kernel, a is the coupling lag or 
phase shift, and the integration is taken over the whole 
volume V. The kernel used is a top-hat kernel with cou¬ 
pling radius R 



Go(r) ^ 



r < R 
r > R 


(3) 


which is normalized as fy Go(r)dr = 1. In simulations, 
the spatial locations are discretized into r = {xi^yi^Zi) 
with 1 < Xi^yi^Zi < L taking integer values in the sys¬ 
tem of linear length L such that V = L^. Hence, the 
control parameters of the system are R and a, with finite 
size effects determined by L. Extensive simulations have 
been done using the Runge-Kutta scheme m with both 
random initial conditions (IC) and specific functions, see 
the Supplemental Material (SM) for details [38]. We use 
periodic boundary conditions (BC) here, yet our findings 
are quite independent of the BC. 

The Kuramoto model with nonlocal coupling is known 
to exhibit chimera states, in which both synchronized 
and unsynchronized regions of oscillators can coexist in 
the same system even though all oscillators are identi¬ 
cal and uniformly coupled. Most studies have been done 
on the one dimensional ring and complex networks |28| . 
In higher dimensions, two qualitatively different chimera 
regimes have been identified for the Kuramoto model 
given by Eqs. 0 , 0 , Eor near global coupling 

with R ^ L and large a < 7r/2, various coherent and 
incoherent strip, spot, plane, cylinder, sphere and cross 
patterns have been observed in two and three dimensions 
(2D and 3D) [39] . The other regime involves shorter 

range nonlocal coupling L R 1 with smaller a. 
In 2D geometries, synchronized spiral waves with unsyn¬ 
chronized chimera cores can appear. They behave like a 
normal spiral, yet the dynamics in the core is unsynchro¬ 
nized [39II42] . Similarly, in 3D, regular scroll waves with 
chimera filaments (or chimera tubes) at their center — 
instead of the linelike filaments of phase singularity — 
have been observed m- 

Existence of knots: Eor L R 1 and large 

effective system size T/R, we observe different stable 
linked and knotted scroll waves in the Kuramoto model as 
shown in Eigs.[^and To clearly visualize the chimera 
tubes and the knotted and linked structures (referred to 
as knots in the following), one has to take into account 
that both phase ^(r, t) and angular frequency C(;(r, t) fluc¬ 
tuate a lot in space as shown in Eig. 03- 0. Thus, 


Figure 1: (color online) Chimera nature of knots {a = 0.8, 
R = 8, L = 200). (a) Snapshot of a Hopf-link shows the 
unsynchronized phase 0(rR) on the isosurface of spatially 
smoothed angular frequency a;(r,t). (b) Mean angular fre¬ 
quency distribution of (ct;(r)), averaged over approximately 10 
periods. The tail to the right corresponds to unsynchronized 
oscillators, (c) Plot of 2D cross-section of ^(r) showing the 
chimera and spiral wave properties. A vertical cut through 
both rings in panel (a), showing four chimera cores, (d) A cut 
through the far edge of a ring, (e) Slicing of a ring, showing 
a ring chimera and two chimera cores of the other ring. 


it is helpful to define a local mean angular frequency 
uj{rR) = fy Go{r—r')uj{r'R)dr. Eig. [la shows a snapshot 
of the chimera tubes by plotting the ^ases of the unsyn¬ 
chronized oscillators for a;(r, t) > const. The presence of 
scroll waves with chimera filaments can also be seen di¬ 
rectly in the phase field. Selected 2D cross-sections of the 
phase field (Eig. Ef) show patterns similar to chimera spi¬ 
rals in 2D uni, while other cross-sections show features 
that are specific to 3D such as the chimera ring shown 
in Eig. [^. Note that the Hopf link and other knots ob¬ 
served generate spherical wave in the far-held. Moreover, 
they are not stationary but keep rotating, drifting, and 
changing their shape over time as shown in the videos 
and the SM [35] . 

As chimera hlaments are associated with scroll waves, 
phase twists can be present along the hlament [43] . This 
is visible in Eig. 01, but can be better visualized by con¬ 
sidering the local mean phase held ^(r, t), dehned by 

= f Go(r - (4) 

Jv 

This is illustrated in Eig. for example. 

Phase diagram: Eor the Kuramoto model given by 
Eqs. Q, Q, and (j^, our numerical simulations allow us 
to obtain a phase diagram as a function of a. This is 
plotted in Eig. [^ together with some of the asymptotic 
states. At small a, only relatively simple scroll wave 
structures with straight chimera tubes are stable. Eor 
ax < a < 0 ^ 0 , also knots such as 1 twist Hopf links 
and 3 twist trefoils (as shown in Eig. [^ are stable over 
hundreds of thousands of scroll wave rotations period T 
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Figure 2: (color online) Phase diagram of the Kuramoto 
model in 3D as a function of a with nonlocal coupling L ^ 
^ 1. Various knots exist between ax and ao. Top panel 
shows stable Hopf-link (left), and trefoil (right) as examples. 
The plots are similar to those in Fig. but smoothed phases 
^(r) are used instead. Shadows on the walls correspond to 
perpendicular projections of the structures. 


123 More stable structures including helices, ring-tubes 
and linked triple rings are shown in the SM [38]. For 
a > 0 ^ 0 , knots as well as simple straight tubes become 
unstable. The evolution in the former case is shown in 
Fig. In the latter case, the dynamics of the chimera 
filament indicates that a finite wavelength instability of 
the filament itself occurs such that the filament grows 
rapidly (see the SM [38]). In both cases, the rapid growth 
of filaments is accompanied by fragmentation through 
collisions leading eventually to an irregular or turbulent- 
like behavior as shown in Fig. Furthermore, in the 
same parameter regime near ao in 2D, chimera spirals 
are stable and no irregular pattern is present [40]. All 
this suggests that the underlying instability is truly 3D in 
nature as the negative line tension instability and similar 
filament instabilities that have been observed in excitable 
and oscillatory media m Fig. also shows that at even 
higher a ^ 7r/2, no filament structures can be recognized. 

The nature of the instability of knots at aK and ao 
are significantly different. Below ax^ any knot trans¬ 
forms through one or multiple reconnections into a single 
untwisted ring which shrinks and disappears, leading to 
homogeneous oscillations. As we have tested, all single 
chimera rings with a radius of up to 80 shrink and eventu¬ 
ally vanish for 0 < a < q^o with no-flux BC (see SM EH]). 
This together with the stable knots for ax < a < ao 
indicates that there is an effective repulsion between fil¬ 
aments in knots that is sufficient to prevent curvature- 
driven shrinkage and stabilize these structures above ax- 
Below ax , the repulsion is too weak to prevent reconnec¬ 
tions. This mechanism is similar to what has been ob¬ 
served for knots in bistable media [12] and plays an im¬ 
portant role in other situations as well El. Simulation 



Figure 3: (color online) Time evolution of Hopf links outside 
their stability regimes near ax and oo, respectively, (a) a < 
ax- After a is decreased slightly from above ax to below, 
the two rings of the Hopf link merge resulting in a single 
ring which eventually vanishes, (b) a> ao. After a is tuned 
slightly from below ao to above, one of the rings grows until it 
collides and reconnects, resulting in a turbulent-like pattern. 


results show that different knots have different stability 
regimes, especially Hopf links are stable over a broader 
range of a than trefoils. Therefore, we denote ax in the 
following as the point at which Hopf links disappear. 

Dependence on R, L, and geometry: Numerical sim¬ 
ulations for 4 < i? < 12 and 64 < L < 300 show that 
the phase diagram presented in Fig. [^is independent of 
the specific choice of R and L as long as L ^ ^ 1. 

Specifically, ax ~ 0.61 and ao ~ 0.90 with uncertainty 
±0.02. The condition L ^ R ensures that finite size ef¬ 
fects do not play a significant role as the size of stable 
knots and the wavelength scale with R [iOllill. For ex¬ 
ample, we find that stable Hopf links cease to exist for 
L/R < 16. Also, if the effective system size is too small, 
more complex knots tend to decay into simpler ones (see 
the SM [^). The condition i? > 1 is also crucial. We 
find that for shorter range couplings R < 3 the lifetime 
of knots is finite [58]. Specifically, no stable knots have 
been observed for local coupling, = 1, independent 
of the IC used to generate Hopf links. This is a con¬ 
sequence of temporal fluctuations in the s hape of the 
individual rings within a Hopf link becoming comparable 
to the minimum separation between the rings such that 
the rings merge and disappear (see the SM [38]) — the 
same behavior as for the instability at ax- We observe 
qualitatively the same for trefoils. 

Robustness with respect to noise: To further quan¬ 
tify the stability of different topological states, we exam¬ 
ine them in noisy environments. This is modeled by an 
additional Gaussian phase noise in the Kuramoto 

model 

0{Y,t) = u}{r,t) + 


(5) 
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Figure 4: (color online) (a) The tolerable phase noise level 
for Hopf links (HL) and trefoils (T). / quantifies the fraction 
of structures persisting after t = 5000 450 scroll wave 

rotations) under noise intensity D. Each point corresponds 
to an ensemble of 10 realizations. Here, a = 0.8, = 4, and 

L = 100. Inset: D* denotes the point of / = 0.5 for a Hopf 
link as a function of R using a — 0.8, L/R — 16. (b) Phase 
diagram showing the stable regime of Hopf links for kernel Gi 
with L ^ R^ 1. The red dot marks the triple point between 
all three phases where ax = <ao- The error bars account for 
different R > 4 and L up to L = 300. Note that 1/ni = 0 
corresponds to the top-hat kernel Go used in Fig. 

where (^(r, t)) = 0 and (^(r, t)^(r', t')) = — — 

As shown in Fig. Hopf links and trefoils can survive 
under noise magnitude as high as 1)* = 0.22. This high 
robustness under noise signifies the topological protection 
of knots. Longer range coupling as quantified by R also 
increases the tolerance of local phase noise as shown in 
Fig. 1^. 

Dependence on spatial kernel: In contrast to the top- 
hat kernel Go, we did not observe stable knots for Gaus¬ 
sian kernels often considered in the context of chimera 
states. This together with the existence of a minimal R 
discussed above indicates that the range of the spatial 
kernel is crucial. To substantiate this further, let us con¬ 
sider the kernel Gi(r) ^ such that Gi ^ Go 

when ni ^ oo\i R < L. Note that ni = 2 is a Gaussian, 
ni = 1 is an exponential, and ni = 0 gives global cou¬ 
pling. Simulations show that if Gi becomes more long- 
ranged as ni decreases, the stable regime ax < < <ao 

of knots shrinks as shown in Fig. I§. This is the only 
effect on the knots as the nature of the associated in¬ 
stabilities along the boundaries appears unchanged and 
follows the pattern shown in Fig. The phenomenon 
is independent of the exact functional form of the kernel 
(see the SM [38]). 

Beyond phase oscillators: The nonlocal CGLE is 

given by m 

A(r, t) = A - (1 + ib)\A\‘^A + (1 + ia)pA, (6) 

where the control parameters are (a, 5). The nonlocal 
coupling Pa is given by 


Figure 5: (color online) (a) CGLE [(a, 6) = (1,0), K = 0.1]: 
Comparison of minimum separation between rings in a Hopf 
link, dmin, and a measure of spontaneous fluctuations in the 
ring shape, (see the SM [38] for details). Bars indicate the 
99% range, (b) Rossler model: Topological superstructure of 
a chimera Hopf link with an attached sychronization defect 
sheet in the period doubled regime (see the SM [38] for more 
viewpoints and dynamics). 

and the coupling strength is K. Since the CGLE can 
be well approximated by the Kuramoto model in the 
weak-coupling limit independent of the specific cou¬ 
pling [25] [45], similar results are expected in certain pa¬ 
rameter regimes. Indeed, stable chimera knots with non¬ 
constant amplitudes |A| exist in the vicinity of the pa¬ 
rameters (a, b) = (1, 0) [46] for K = 0.1 (see the SM [^ ) 
and larger values of K. The lifetimes of knots are longer 
than 6 x 10^, provided that R^ 1. All results discussed 
above for the Kuramoto model also hold qualitatively for 
the CGLE. This includes in particular the break-up of 
knots for small R. Fig. provides a clear rationale why 
this happens: The separation between the rings in a Hopf 
link shrinks with decreasing R such that it eventually be¬ 
comes comparable to the amplitude associated with the 
temporal fluctuations in the shape of individual rings. 
This offers an explanation of why no stable knots have 
been observed in the CGLE with local coupling. Our 
findings for the CGLE imply that all oscillatory systems 
with appropriate nonlocal coupling should exhibit stable 
knots in some parameter regime near their Hopf bifurca¬ 
tion. 

Gomplex oscillatory systems: We also observe stable 
chimera knots if the uncoupled oscillators are far from the 
Hopf bifurcation and undergo complex or even chaotic 
oscillations, requiring at least a three-dimensional local 
phase space. A specific exampie is the Rossier modei with 
noniocai coupiing 1171 , which exhibits a phenomenoiogy 
with many features in common with those observed in 
compiex osciiiatory systems inciuding chemicai experi¬ 
ments [48]. It is given by 

X(r,t) = -Y-Z+px, 

y(r,t) = X + aF+py, (8) 

Z(r,t) = b + Z{X-c), 


PA{r,t)=K j Go(r-r')[A(r',t) - A(r,i)]dr', (7) 


where the controi parameters are (a, 5, c) and the nonio¬ 
cai coupiing px(r, t) and py (r, t) are defined anaiogousiy 
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to Eq. 0. For a = 6 = 0.2, the effective a decreases as 
c increases For R 1, we observe stable chimera 

knots in the period-doubled regime (c = 3.6) and in the 
chaotic regime (c = 4.8) with weak coupling K = 0.05 
(see the SM [38]). All findings described above for the 
other models hold qualitatively here as well. Stable knots 
only exist if the coupling between the oscillators is nei¬ 
ther too short-ranged nor too long-ranged. For exam¬ 
ple, we did not observe stable Hopf links or trefoils for 
R = 1 OT when the kernels were Gaussian in the parame¬ 
ter regimes given above. Moreover, in the period-doubled 
regime, synchronization defect sheets (SDS) — the ana¬ 
log of synchronization defect lines in 2D systems [47l[49| 
— can be observed for the first time and, more impor¬ 
tantly, connect the different filaments (see Fig. [^, the 
SM and movie [38]). This leads to another layer of topo¬ 
logical structure associated with the knots, making this a 
unique phenomenon and adding potentially to their gen¬ 
eral robustness if multiple knots are present |48) . 

Discussion and conclusions: Our findings show that 
knots exist and are stable over a significant range of pa¬ 
rameters in various oscillatory systems with nonlocal cou¬ 
pling as long as the characteristic coupling length of the 
kernel is sufficiently large and the tail of the kernel decays 
sufficiently fast. The variety of knots is also much higher 
compared to what has been reported for excitable me¬ 
dia miiia. For example, we have also observed other rel¬ 
evant unknotted structures such as stable double helices 
(see the SM [38]) — a structure that has remained elusive 
in the CGLE with local coupling m- This suggests that 
the models considered here can serve as paradigmatic 
models to study various knotted and unknotted struc¬ 
tures associated with scroll waves in general, including 
the novel topological superstructures of knots with SDSs. 
More specifically, it allows one to explore the topological 


constraints imposed by the phase field on the observable 
phase twists associated with a given knot — a field largely 
untouched [33] — as well as the effect of synchronization 
defect sheets on knots for the first time. 

A remaining open question is to which extent the ex¬ 
istence of stable knots in oscillatory systems depends on 
the presence of a chimera state. While our findings sug¬ 
gest that a chimera state is a necessary condition, there 
is no fundamental reason to substantiate this. However, 
our simulations indicate that the mobility of the scroll 
wave filaments plays an important role. If the filaments 
move or meander sufficiently fast (e.g. = 1 or for a 

Gaussian kernel with large a), no chimera state can be 
numerically observed and stable knots are absent. This is 
similar to what has been reported for chimera spirals in 
2D [41] and knots in excitable media [TT]. One possible 
way forward is the recently proposed ansatz by Ott and 
Antonsen ED ED, which has been successfully applied to 
study the existence and stability of chimera spirals |53| . 

In addition to the robustness of knots under dynami¬ 
cal noise, we also find that Hopf links and trefoils can 
emerge in a self-organized way from random IC with 
fair probability (see the SM [38]). Both features indi¬ 
cate knots should be observable in real-world oscillatory 
systems that follow a dynamics similar to the models 
studied here, with most likely candidates to be chemi¬ 
cal systems 123 mi ED ED Eg . Yet, the observation of 
chimera filaments in natural systems remains a challenge 
for the future. 
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Supplemental Materials: Linked and knotted chimera filaments in oscillatory systems 


TOPOLOGICAL STRUCTURES 


Fig. 


ST] shows various long-lived (stable or metastable) topological structures in the non-local Kuramoto model 

Note that knotted structures more complicated than simple Hopf links tend to have 


within the regime ax < a < ao. 
smaller stability regimes. An exception are (knotted) structures that require periodic boundary conditions (BC) and 
do not drift, which can be stable below ax- This is shown, for example, in Fig. [^and includes straight filaments. 
Simulations suggest that the multi-filament structure in Fig. [S^ is stable for a > 0. 


II. DYNAMICS 


In the non-local Kuramoto model, knotted structures that exist independent of the specific choice of BC (periodic 
vs. no-flux) are not stationary but drift, rotate and change their shape over time. As an example. Figs. [S^a) and 
[S^b) show a few snapshots for different structures. Note that for the system sizes studied, the center of mass motion 
is not straight over long time scales. The phase field away from these knotted structures takes on the form of spherical 
waves as shown in Fig. S4 In case of the ring-tube structure (which is specific to periodic BC), the ring propagates 
along the tube and keeps distorting the local part of the tube while it travels, see Fig. S3 (c). In all these cases, the 
direction of the filament motion can be deduced from the instantaneous angular frequency uj{x^y^z) shown in the 
rightmost column of Fig. [S3| 


III. CREATING CHIMERA KNOTS 
A. Random initial condition 

Knots and links can appear spontaneously from random initial conditions (IC). The transient time is of the order of 
one thousand scroll wave periods in the regimes being studied. A few snapshots of typical transient states are shown 
in Fig. |S5[ Using random IC, we can obtain all knotted structures shown in Figs. |S1^-|S1|F. The specific probabilities 
of generating Hopf links and trefoils from random IC are summarized in Table 


B. Algorithm to create rings and Hopf links 

First, the phase field of a single ring is considered. Suppose the center of a ring is located at tq = {xo^yo, zq) with 
radius Rq and the normal vector of the ring is pointing in the positive z direction. A parameterization of the location 
of this ring using (j) G [0, 27r) is 


r(0) = {rx,ry,rz) = {xo + Rocos0,^o + Rosmcf), zq). (SI) 

To create a ring shaped filament corresponding to phase singularities, the phase field needs to be specified in the 
whole domain such that it is smooth outside the ring but results in 27r phase difference while going around a point 
on the filament. This can be done by defining 

(^) 

/ = y(a; - + (y - 2/o)^ (S3) 

for any spatial point r = (x,^, z). Then the phase of each oscillator 0{t) can be computed by 0{t) = '0(r), where 

'0(r) = kd — if — s(j) — (S4) 


Here, k is the wavenumber, d is the distance to the closest point on the ring d = {x — + (^ — TyY ^ {z — 

ip is the angle between the plane consisting of the ring and the line to the closest point of the ring, s is the twisting 
number, (p is the ring parameterization, and /3 is a constant phase shift. Examples are shown in Fig. |S6[ 
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Figure SI: Non-trivial topological structures in the non-local Kuramoto model with periodic BC for L = 100. The lifetime 
T and the corresponding parameters are given for each subfigure, r > tq means that the structure is stable within the testing 
time limit tq, while r ~ tq means the structure breaks down around tq (order of magnitude). The period of the scroll waves is 
about T ~ 11 for a = 0.8. This implies a lifetime of more than lO^T for Hopf links and trefoils. Together with the robustness 
in the presence of noise as established in the main text, this suggests that the lifetime r ^ oo when L ^ R. 
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The phase field of a Hopf link can be created by combining two rings, requiring a method to smoothly superimpose 
them. This can be achieved using a distance dependent phase: 

^(r,ro,s,^)= (S5) 

which is based on the inverse square distance. Then the phase field of a Hopf link 0{t) can be calculated by 

p(i.)e*®W = ^{x, y, z, xq - i?o/2, yo, zo, s,/3 = 0)+ ^(x, z, y, xq + Ro/2, yo, zq, s, ^ = tt) (S6) 

with twisting number 5 = 1. Examples are shown in Fig. |S7| Note that a structure in a given system size L can be 
rescaled to L' using a simple scaling function of the form 9'{x'^y'^ z') = , [^y'\ , where [-J denotes 

the floor of the number (which is necessary since the oscillators are arranged on a discrete lattice) and the prime 
denotes the new phase and new location. This rescaling works quite well for the top-hat kernel as long as i? ^ i?' ^ 1. 
Also, if the smoothed phase 0 of knots — see Eq. (4) in the main text — is used as IC, the unsynchronized region 
around the filaments can redevelop. 


C. Reconnecting chimera filaments using random patches 

A new structure can be obtained by reconnecting local filaments of a known structure. This reconnection requires a 
detailed specification of the whole local phase field that is smooth, without creating new filaments and while matching 
the desired filaments. This can be hard to do if the local filaments are obtained from a simulation. Alternatively, 
based on the observation that only simple straight filament can form in a small system size L/R from random IC, it 
suggests a way to transform a structure by randomizing a whole local region. Using this method, we have successfully 
created trefoils and a few other knots. To begin with, a structure that is similar to the desired knot is needed, with 
the region of reconnection close to each other. For example, the structure in the red box shown in Fig. is a trefoil if 
the top parts are connected. After half a dozen trials using different shapes of the randomized region, we were indeed 
able to create a trefoil. Note that the region should be large enough to form a tube but not too large to form other 
structures. This method may suggest a similar way to create knots in real world experiment. 


IV. INSTABILITIES 

In the main text, the instabilities at ax and ao of Hopf links in the Kuramoto model have been discussed. The 
instability near ax is caused by a lack of repulsion to counter curvature-driven shrinkage, so knots collapse and 
disappear. On the other hand, the instability near ao originates from an instability of the filament where the filaments 
become longer and longer and eventually collide with themselves or other filaments. This effect is particularly clear in 
large domains as shown, for example, in Fig. Note that an elongation also happens as a transient state when the 
parameters are suddenly changed or starting from a non-perfect IC. However, it will eventually shorten after refolding 
to an asymptotic state as also observed for other models |S1| . Other instabilities are discussed below. 


A. Instability of a single ring 

Direct simulations show that rings are not stable for a < ao with no-flux BC. As shown in Fig. all rings shrink 
in size and eventually vanish. The largest ring tested had radius Ro = 80. This shrinkage process occurs for both 
nearest neighbor coupling R = 1 and non-local coupling R = 4. Note that the time it takes for a ring to disappear is 
approximately the same in both cases for the same effective radius Ro/R and effective system size L/R. Also, almost 
all transient (knotted) states resulting eventually in homogeneous oscillations become rings in their penultimate stage. 


B. Instability of knots for R—1 


As shown in Fig. |S7[ using the IC for Hopf links described in Section HIB| can result in a stable knot if R ^ 1. For 
the choice of Rq, the two rings initially shrink in size and then an effective repulsion prevents further shrinkage. At 
the same time, the center of the Hopf link starts moving. In contrast, if the same IC is used with nearest neighbor 
coupling R = 1, the two rings will eventually collide with each other and decay into a single ring, which in turn shrinks 
and vanishes. 
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C. Filament instability at 

As illustrated in Fig. |S9[ the instability at o^o for simple straight chimera filaments is characterized by the emergence 
of secondary structures and the elongation of filaments. The same qualitative behavior is observed for knotted 
structures in the regime a > q^o as shown in the main text. Nevertheless, the knotted structures can persist for 
thousands of scroll wave periods before they break up consistent with critical slowing down near a phase transition. 


D. Instabilities from finite size effects 

A stable knotted structure becomes unstable when it is confined in a small effective system L/R. While a Hopf 
link simply decays into a single ring which eventually vanishes, the situation is more complicated for larger and more 
complex knotted structures. One example is shown in Fig. |S10| starting from a triple ring for = 4 in L = 100, which 
decays into a ring knotted with 8-shape ring, and then transforms into a trefoil. Depending on the IC and the exact 
parameter regime, the decay path can be different. Note that even for the moderately larger system size L = 150, 
triple rings have significantly longer lifetimes (r > 20000) in some parameter regimes. 


V. SPATIAL KERNELS 


As mentioned in the main text, our main findings do not depend qualitatively on the exact functional form of the 
considered kernels. For example, using the kernel 


1, \x\,\y\,\z\ < R 

0, otherwise 


(S7) 


instead of the top-hat kernel Go gives pretty much identical results for the stability of knots. As another example, 
using the kernel 


G2(r)-(l + e"^(’-^))-i (S8) 

with an exponential tail instead of the kernel Gi with super-exponential tail exhibits the same phenomenology: With 
decreasing n 2 , the stable regime of knots shrinks. The shape of the kernels at the transition points of Gi and G 2 for 
R = A and a = 0.8 are shown in Fig. |S11[ 


VI. OTHER OSCILLATORY MODELS 
A. Non-Local Complex Ginzburg-Landau equation (CGLE) 

The non-local CGLE considered here is [S2]: 


A(r, t) = A — (1 + ^6)1^1^^ + 


K(1 


«)/ 


G(r-r')(A(r')-A(r))dr', 


(S9) 


where the control parameters are (a, 6), the coupling strength is K and G = Go in the following. Under sufficiently 
weak coupling K ^ 0, the local field oscillates with unit amplitude |A| ^ 1 and behaves like a simple phase oscillator 
in the non-local Kuramoto model. Therefore, we can use the knotted structures found in the Kuramoto model as IC 
by simply setting A(r, t = 0) = We find that one of the r egimes with stable Hopf links is 0.95 < a < 1.15 

for 6 = 0 and K = 0.1 provided that L ^ R ^ 1. In Fig. S12[ a), the phase portrait shows that the magnitude 


of all oscillators onl y dev iates slightly from \A\ = 1 in this case. For stronger coupling K = 0.2, the deviations in 
A increase (see Fig. |S12| ^b)) but stable knots still exist. In bot h ca ses, the pha se 0{r) = arg(A(r)) behaves similar 
to the Kuramoto model as confirmed by Fig. |S13[ d). As Figs. |S13[ b) and |S13[ c) show, the chimera nature is also 
evident from the Ke{A{x^y,z)) and lm{A{x^y,z)) fields. Using the local mean field ^(r), one can easily locate the 
unsynchronized filaments |S6| . An example is shown in Fig. |S13| ^a). 
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B. CGLE: Minimum separation &; spontaneous fluctuations 


When R becomes too small, knots are no longer stable. This instability can be characterized by the dynamics of 
the filament(s) that make up the knots. Even though the region around the filament is unsychronized, the filaments 
can be found by a filament detection algorithm |S3| of the mean field (see Fig. |S14[ a)). The length of filament can 
therefore be defined as the number of occupied lattice sites. Denote the two r ings or filaments of a Hopf link as 
Fi and F 2 with circumference (or length) Ci and C 2 , respectively. As Fig. |S14[ c) shows, Ci and C 2 fluctuate over 
time in a synchronous way. Fluctuations are also present in the minimum separation between Fi and F 2 , defined 
as dmin = minriGFi ( 1 * 1 , r 2 ), as shown in Fig. |S14[ e). To characterize these fluctuations statistically and identify an 
associated length scale, we consider the difference between the 99.5%-quantile and the 0.5%-quantile associated with 
dmin^ corresponding to the error bars shown in Fig. 5(a) in the main paper. As shown in Fig. |S15[ b), this difference 
is not varying much across the considered values of R. This is in sharp contrast to the linear scaling of d^in with R 
(see Fig. 5(a) in the main paper). 

To substantiate that the intrinsic length scales associated with filamen t fluc tuations do not strongly vary with i?, 
we further consider the fluctuations of a single straight filament (see Fig. |S14| ^b)). To ensure a fair comparison with 
the fluctuations of Hopf links, we choose a system size L = Tg// such that the average single filam ent length {£) equals 
the average circumference {C) = {{Ci) + (C 2))/2 of the fil amen ts in the Hopf link (see Fig. S14[ d)). The dependence 
of both these quantities as a function of R is shown in Fig. |S15[ a). To characterize the fluctuations of a single straight 
filament, we calculate its roughness. Due to the chosen initial conditions, the roughness is identical to the deviation 
from a straight filament oriented along the ^-axis. Specifically, we define the deviation from the straight filament 
center r^y = (1/L) ^xy{z) to be 




'■xy 


(^) 


'■xy\^ 


(SIO) 


where Txy{z) is the intersection point of the filament with the x-y plane for a given 2 ;. The roughness /S^xy = 


{1/L)^z ^xy{z) is now simply Axy{z) averaged over z. As Fig. S14T) shows, the roughness varies over time. To 


characterize these (non-negative) fluctuations in the roughness over time and within an ensemble and to identify an 
associated length scale, we consider the 99%-quantile and denote it by This is the quantity shown in Fig. 5(a) in 
the main paper and again it does not vary much across the considere d va lues of R. For a direct comparison with the 
length scale of fluctuations in the case of a Hopf link, please see Fig. |S15[ b). 


C. Non-Local Rossler model 

The non-local Rossler model considered here is [S4] : 


X{v,t) 

= -Y-Z+KJ 

' G{r-r'){X{v')-X{v))dY, 

(Sll) 

Y{T,t) 

= X + aY + K J 

''G{r-r'){Y{r')-Y{r)) dr', 

(S12) 

Z{v,t) 

= b + Z{X-c), 


(S13) 


where the control parameters are (a, 6, c) and the coupling strength is K. Again, we can use (X, T, Z) = (cos 6>, sin6>, 0) 
with 0{t) from states with knotted structures generated by the Kuramoto model as IC. When a = 6 = 0.2, the effective 
\a\ decreases as c increases |S5| . We observe stable knots within 3.3 < c < 5 for weak coupling K = 0.05 provided that 
L ^ R ^ 1. Note that as c increases the intrinsic dynamics of the oscillators also changes. Namely, the dynamics 
undergoes a period-doub ling c ascade to chaotic oscillations. In particular, we o bserv e stable knots in the period-2 
regime with c = 3.6 (Fig. |S1^ , as well as in the chaotic regime with c = 4.8 (Fig. |S17[) . 

An additional feature of the wave dynamics in these regimes is evident from Fig. |S16| and Fig. |S17( The amplitudes 
are modulated. For example, in the period-2 regime alternating wave maxima are present. A topological consequence 
of such a behavior is that two dimensional structures exist such that the local dynamics has a lower period than 
that of the bulk. Specifically, these structures, called synchronization defect sheets (SDSs) in the following, separate 
domains of different oscillation phases and for periodic BC either originate from a filament or are closed. More 
importantly, every filament has an attached SDS such that they become part of any knotted structure. This can 
already be observed in the cross-sections shown in Fig. |S16| and Fig . |S17| To clearly identify SDSs, we use the 
detection algorithm developed for the lower dimensional case in Ref. |S5|. A specific example of SDSs is shown in 
Fig. |S18| and their subsequent motion is shown in the Supplementary Video. 
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Figure S2: Various long-lived filaments with periodic BC, L — 100 and = 4. (a-b) Each filament connects with itself through 
one of the surfaces, (c) Each filament passes through all three surfaces before connecting back to itself. 


Table I: Spontaneous formation of Hopf links and trefoils in simulations with random IC and periodic BC. 


a 

R 

L 

Number of simulations 

Number of Hopf links 

Number of trefoils 

0.7 

4 

100 

200 

1 

1 

0.7 

5 

100 

500 

4 

0 

0.8 

4 

100 

500 

13 

0 

0.8 

5 

100 

500 

5 

0 

0.7 

5 

200 

100 

3 

1 

0.8 

8 

200 

100 

4 

0 


























t = 0.0 t = 50.0 t = 400.0 



(a) Hopf link 



(c) Ring-tube 
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100 




Figure S3: Dynamics of different topological structures {a = 0.8, = 4 and periodic BC). The first three columns are snapshots 

at three different instances in time. The rightmost column is the iso-surface plot of the instantaneous angular frequency u of 
the last snapshot. Blue indicates the region with |a;| < |a;|, while orange indicates the region |ct;| > |ct;|. Filaments are moving 
away from the orange region. 














































9 



Figure S4: The spherical wave generated around a Hopf link is shown by plotting the iso-surface 0 — ^ {a — 0.8, R — S and 
periodic BC). The irregular pattern at the center is the region with unsynchronized oscillators that form the filament. Note 
that only the lower half of the system is shown to highlight the structures near the center. 



200 

(a) A transient state for R = S, L = 200. 


0 



(b) A transient state for R = A, L = 200. See 
Sec. |III C] for a discussion of the significance of 
the red box. 



Figure S5: Some snapshots of transient states {a = 0.8 and periodic BC). 
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Figure S6: Shrinking rings with (a) non-local coupling = 4, (b) nearest-neighbor coupling R = 1. Parameters: LjR — 50, 
R^ijR — 20, a = 0.8 with no-flux BC. 
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Figure S7: Formation of knots for (a) non-local coupling R = 4, (b) nearest neighbor coupling R = 1. Note that the IC are 
exactly the same in both cases. Parameters: Ro = 25, L = 100, a = 0.8 with no-flux BC. 



Figure S8: Instability of a trefoil near the transition point ao (i? = 4, L = 150 and periodic BC). This snapshot shows the 
initial elongation of one branch of the trefoil, which has collided with itself and formed an extra ring. 
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0 0 0 



Figure S9: Snapshot series of the instability of straight filaments at a = 0.95 > a^. (a) Single filament with L = 100, = 4 

and no-flux boundary conditions. Some secondary structures develop with local twisting before break-up. (b) Two filaments in 
a larger domain L = 200 and R — A with PBC. The rapid elongation of one of the filaments is evident. 



Figure SIO: Decay of a triple ring for L — 100, R — a — 0.7 with PBC. (a) t = 0, triple rings, (b) t — 15000, decay into a 
ring knotted with an 8-shape ring, (d) t = 20000, further decay into a trefoil. 
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Figure Sll: (color online) Plot of the localized kernels that can still result in a stable Hopf link, with estimated critical values 
ni = 3.8 and n 2 = 1.9 for = 4 and a = 0.8 of Kuramoto model. Go is shown for comparison. 




- 1.0 - 0.5 0.0 0.5 1.0 

Re(A) 


(a) K = 0.1 


(b) K = 0.2 


Figure S12: Snapshot of the states of the oscillators in phase space for a Hopf link in the non-local CGLE for (a, 6) = (1,0), 
L = 200, R = S and periodic BC. 
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(b) Re(A) (c) Im(A) (d) 6 = arg(A) 


(a) shows the unsynchronized region 
An x-y cross-section of the different fields at z = 100 is plotted in (b)-(d). In (b) and (c), 


Figure S13: Snapshot of a Hopf link in the non-local CGLE corresponding to Fig. |S12 
corresponding to the chimera knot 


the color map from deep blue to red corresponds to values from —1 to 1 in the respective field. 
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(a) Snapshot of the filaments of a Hopf link. 


0 



(b) Snapshot of a straight filament with no-flux BC. 



(c) Circumferences Ci of the two rings shown in (a) as a (d) Length of the filament i shown in (b) as a function of time, 

function of time. 



(e) Minimum separation between the two rings shown in (a) as 
a function of time. 



t 

(f) Roughness of the filament shown in (b) as a function of time. 


Figure S14: CGLE with = 4, a = 1, 5 = 0, iC = 0.1. (left) Temporal evolution of a Hopf link with L = 80. (right) 
Temporal evolution of a single filament oriented along the z-direction with L = 91 which results in a time average filament 
length {£) 110 that is approximately the same as the time average circumference (C) ~ 110 of the rings in the left column. 
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(a) Average circumference and effective system size L^ff. 



(b) Comparing different length scale of fluctuation. 


Figure S15: CGLE with a = 1, h = 0, iC = 0.1 as in Fig. |S14| (a) System size, Fe//, for which a single filament has the 
same average length {1) as the average circumference {C) of a Hopf link, (b) Measures of fluct uations for the case of a Hopf 
link (spread in the minimum separation betw een t he rings, dmin,99.5% ~ dmin,0.5%^ see Fig. S14^e)) and for the case of a single 
filament (spread in the roughness, see Fig. |S14[f)), both as a function of R. 
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(d) Z (e) 6> = tan-i(y/X) 


Figure S16: Snapshot of a Hopf link in the non-local Rossler model for (a,6, c, iF) = (0.2,0.2,3.6,0.05), corresponding to the 
period-2 regime. The lifetime of the knot is r > 10^. The same 2D cross-sections of the Hopf link are shown in (b-e) for the 
different fields X^Y^Z and 0. The color scheme is such that deep blue represents the most negative value, and red represents 
the most positive value. The discontinuities of color along the wave fronts correspond to cross-sections of synchronization defect 
sheets. L — 200, R — S and periodic BC. 
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(d) Z (e) 6> = tan-i(y/X) 


Figure SI7: 


Similar to Fig. S16 


but in the chaotic regime (a,5, c, iF) = (0.2,0.2,4.8,0.05). 
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(f) z = 109 (g) z = 118 (h) z = 130 

Figure S18: Visualization of the synchronization defect sheets (SDSs) present in Fig. |S16[ (a), (b): Different 3D plots of the 
SDSs. (c)-(h): 2D cross-sections at different values of z. The red lines represent the cross-sections of SDSs and the yellow-red 
dots indicate the unsynchronized regions. Note that the cross-section in (e) is the same cross-section as in Fig. |S16| (b-e). 














